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ABSTRACT 

The morphology of galaxy clusters reflects the epoch at which they formed and hence 
depends on the value of the mean cosmological density, JIq. Recent studies have shown 
that the distribution of dark matter in clusters can be mapped from analysis of the 
small distortions in the shapes of background galaxies induced by weak gravitational 
lensing in the cluster potential. We construct new statistics to quantify the morphology 
of clusters which are insensitive to limitations in the mass reconstruction procedure. By 
simulating weak gravitational lensing in artificial clusters grown in numerical simula- 
tions of the formation of clusters in three different cosmologies, we obtain distributions 
of a quadrupole statistic which measures global deviations from spherical symmetry 
in a cluster. These distributions are very sensitive to the value of fig and, as a result, 
lensing observations of a small number of clusters should be sufficient to place broad 
constraints on VLq and certainly to distinguish between the extreme values of 0.2 and 1. 
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1 INTRODUCTION 



Forman 199E . Motivated by these interesting but contro- 



versial observ ations suggesting recent c luster growth, 



In a lo^' density universe, density fluctuations cease to grow et al. (1994) and Mohr et al. (1995) demonstrated that 



after a redshitt z ~ 1, 



where S 2o is the p resent value of 
[Peebles 



Evrard 



1980 



the cos nological density parameter (see e.g. 
§11 & §13). Introducing a cosmological constant, Ap (which 



in agreement with analytic predic tions ( Richstone, Loeb, 
[Lacey fc Cole 1993D , 



Turner 1992 



the internal structure of 



cluste rs is sensitive to the cosmological model. 



Evrard et al 



we sha jl express in units of SHq, where Hp is the present (1994| ) constructed X-ray surface brightness maps from their 

hydrodynamic N-body simulations of cluster formation and 
compared these to X-ray maps of 65 clusters observed with 
the Einstein Imaging Proportional Counter. They concluded 
that galaxy clusters with the observed range of X-ray mor- 
phologies were more likely to have arisen in a high f2o cos- 
mology. 

In parallel with these developments, it has become clear 
that the structure of the dynamically dominant dark matter 
in clusters can be reliably mapped by analysing the distor- 
tions in the images of background galaxies induced by weak 
gravitational lensing in the cluster potential. Our aim in 
this paper is to assess whether the surface overdensity maps 
revealed by such analyses may be used to quantify cluster 
morphology and hence to constrain f2o. To this end we sim- 
ulate the lensing of background galaxies by clusters drawn 
from three different cosmological mode ls. These are the same 
clusters studied by Evrard et al. (1994 ). Fo r each cluster we 
create a mock CCD frame as described in Wilson, Cole, & 



value of the Hubble constant) makes this cessation more 
abrupt. Hence, in low fio models, both with and without 
a cosmological constant, clusters form at moderate redshift 
[z = 1-4), and subsequently accrete very little material. In 
this period, which can span many dynamical times, the in- 
ternal structure of the clusters relaxes to produce smooth, 
nearly spherical configurations. In contrast, in Slo = 1 mod- 
els structure formation occurs continuously, rich galaxy clus- 
ters form in abundance only at very low redshift [z — 0.2- 
0.3), and continue to accrete material even at the present 
epoch. Hence, if Qo ~ 1, many clusters are expected to show 
evidence of recent merger events and to have irregular mor- 
phologies. 

Over the past few years a number of optical and X-ray 
studies have suggested that a significant proportion of clus- 
ters , _2erha£s_j;;;^^0% , slww_eyid^ 



Geller fc Beers 1982|; bressler fc Schectman 19881 IWest & 



Bothun 1990t [Forman fc Jones 1990| ; [Bird 1994 [West, Jones 
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Frenk (1996) (hereafter WCF) and then use the inversion 



technique of Kaiser & Squires (1993) to reconstruct a map 



of the projected cluster overdensity. 

The remainder of this paper is organised as follows. 
The main features of the Kaiser & Squ ires reco nstr uction 
method are summarized in Section 2.1, Section |2.2 intro- 



duces new statistics for measuring cluster shapes, specifi- 
cally tailored to be insensitive to uncertainties inherent in 



2.3 



the reconstructed surface overdensity maps. In Section 
we describe the simulated galaxy clusters to be used as lenses 
and whose density profiles we examine. Section reviews 
how we generate the background distribution of galaxies and 
create mock CCD frames of the lensed galaxies. Our results 
are presented in Section^ where we compare the reconstruc- 
tions to the original clusters and assess how the measured 
distribution of cluster shapes depends on the cosmological 
model. We conclude in Section ^ with a discussion and sum- 
mary of our main results. 



2 METHODS 

2.1 Lensing Reconstruction 

The lensing reconstruction technique employed in this paper 



was devised by Kaiser & Squires (1993) (hereafter KS). It 



produces a map of the estimated overdensity, o-{d), where 
a{d) is the deviation of the true cluster surface density, S, 
from the mean surface density, S, within the area being 
considered, measured in units of the critical surface density 



Sc: 



s-s 

'S'crit 



(2.1) 



where is the position vector on the plane of the lens. 5crit 
is defined as 



•S'crit 



Do 



4tvG DoiDi, 



(2.2) 



where D denotes angular-diameter distance and the sub- 
scripts refer to the observer, lens and source galaxy. This 
critical surface density is the minimum required to produce 
multiple images of a source object. 

As shown in WCF, the KS technique recovers the mass 
distribution in clusters extremely well when the lensing is 
weak i.e. when the bending angle varies slowly with po- 
sition. Our CCD simulations in that paper demonstrated 
that observational complications such as seeing and noise act 
to diminish the reconstructed surface density. However, we 
showed that it was possible to correct fully for this diminu- 
tion and proposed a method to estimate a multiplicative 
compensation factor, f. 



2.2 Quantifying Shapes 

Lensing reconstructions such as that of KS reliably return 
relative values of the surface overdensity, a{6) (equation 2.1), 
but do not fix its absolute value. Since 5crit depends on the 



geometry of the lensing configuration through the angular- 
diameter distance relationship, the redshift distribution of 



source galaxies must be known before the mean value of 
S'crit (equation 2.10 below) and hence the absolute surf ace 

In 



2.1 



overdensity S — S can be obtained from equation 
addition, the KS technique is only sensitive to variations in 
surface density from the mean. This is because a uniform 
slab of material across the whole lens plane does not distort 
the images of galaxies lying behind. Thus, the mean surface 
density, S, is also unknown unless the region analysed is suf- 
ficiently large to encompass the whole of the lensing material 
so that the surface density near the edge of the region can 
be taken as the zero-point. 

In order to minimise the effects of these uncertainties, 
we define dipole and quadrupole statistics which depend not 
on the absolute surface overdensity, but only on the values 
of d-{9) relative to each other. These statistics, therefore, 
are explicitly independent of S, S'crit and any compensation 
factor /. We define the dipole, D{A), by 



D{A) = 



[dj + dl 



where 



d2 



_\di 

Ai 

)xdA, 
)ydA 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



Here, tTcon is the surface overdensity of the contour level 
above which we choose to evaluate the statistic and A is the 
area within that contour. H is the Heaviside step function, 
H{cT — aeon), which is equal to if a < (Jcon and equal to 1 
if (T > (Tcon • The integrals are evaluated over the area within 
the contour. Similarly the quadrupole, Q{A), for the same 
area A is defined as 



and 

A^ I H{a - acon)dA. 



where 



A^ 



M 

A^ 



qi ^ I H [a ~ Gcon){x^ - y^)dA 



and 

52 = / H{a - acon)'2xydA. 



(2.7) 



(2.8) 



(2.9) 



The X and y coordinates are measured relative to the cluster 
centre. (We shall discuss exactly how we define the cluster 
centre in Section |^.) In practice, the integrals become sums 
over pixels in the reconstructed maps. Note that both these 
statistics are dimensionless. If one pictures a reconstructed 
surface density map as a contour plot, then these statistics 
are independent of the value of a labelling each contour 
because we select the contour by the area A that it encloses 
rather than by its level. Thus, once a suitable area A has 
been chosen, say 5 arcmin^, the statistics are independent 
of :- 

(i) Any dilution of the lensing signal due to seeing. 
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(ii) Nonlinearity effects which may result in an underes- 
timate of the surface overdensity in regions of high 5*. 

(iii) The (unknown) value of Scrit. 

The values of the statistics D{A) and Q{A) will, of course, 
depend on the noise level in the reconstructed surface den- 
sity maps. It is this effect on these statistics that we aim to 
quantify with our simulations of lensing. 



2.3 The Cluster Simulations 

As our lenses we used a set of eight N-body gasdynamic sim- 
ulations of the formation of galaxy clusters, further details 



of which may be found in Evrard et al. (1994). The clusters 



evolved from the same eight sets of initial density fields but 
in three different cosmologies: 

(i) A biased Einstein-de Sitter model [ f2o = 1, ""8 = 0.59, 
where as is the rms fluctuation of mass within spheres of 
radius 8/i~^Mpc, and h is the present value of the Hubble 
constant in units of 100 kms~^Mpc~^]. 

(ii) An unbiased {ag = 1.0) open model with fio = 0.2 
and Ao = 0. 

(iii) An unbiased (as = 1.0) low density flat model with 
Qo = 0.2 and Ao = 0.8. 

In each of fou r periodic box es of size L = 15, 20, 25 

created two constrained 



Evrard et al. 



and 30 h'^ Mpc, 
realizations of Gaussian random density fields, with the 
standard cold dark matter power spectrum appropriate for 
Ho — 50 km s""*^ Mpc ""*^. In each case, using the technique 



of Bertschinger (1987), they imposed the constraint that, 
when smoothed with a Gaussian of width 7?/ = 0.21/, there 
be a peak at the centre of the box with height 2.5-5 times 
the rms density fluctuation on this scale. 

Evrard's (1988) P^M/SPH hydrodynamic N-body code 
was used to evolve the particle distributions to the present 
epoch. Two sets of 32^ particles represented the dark mat- 
ter and gas respectively. A baryon content of Sib =0.1 was 
assumed for all the models. Gravity, PdV work and shock 
heating were incorporated for the gas but the effects of ra- 
diative cooling were ignored. The spatial resolution was ap- 
proximately 0.005L, varying from 75 to 150 kpc depend- 
ing on the box length. For our lensing simulations we chose 
the output epoch from each cosmology closest to redshift 
z — 0.18. This value is fairly typical in observational studies 
(Wilson et al. 1996) and corresponds to the value adopted 
in our previous simulation s (WCF). 



In the simulations of Evrard et al. (1994) the mass of 



each cluster was approximately proportional to the volume 
of the simulation box. Our aim here is to study the morphol- 
ogy of clusters detected by gravitational lensing with similar 
signal-to-noise ratio in each cosmology. Thus, for our pur- 
poses it is more convenient to have a set of clusters all of the 
same mass. To achieve this we simply rescaled the mass of 
each particle in the smaller clusters so that the total mass 
was the same as that in the largest cluster i.e. the one in 
the L — 30h~^ Mpc box. This required multiplying each 
mass by the factor (30/i~^Mpc/L)'^. In order to preserve the 
correct density, we also multiplied the simulation length by 
30/i"^Mpc/L. 



The typical mass of a rich cluster is known approx- 
imately from observations. For a fair intercomparison of 
the models, the clusters considered should all have approx- 
imately the same mass within a flducial radius such as the 
Abell radius (1.5/i~^Mpc). The simulated clusters in the 
open and flat Qq = 0.2 cosmologies were somewhat less 
massive and so we scaled them as in the previous paragraph 
to have approximately the same total mass as the average 
fio = 1 cluster (~ 1.36 x lO^^/i'^M© within an AbeU ra- 
dius). The maximum scaling required was a factor of 2.4 in 
mass. 

After this rescaling the initial conditions are no longer 
appropriate for a standard CDM model but instead have 
power spectra with slopes which differ slightly from the cor- 
rect ones. However, analytic arguments show that moderate 
changes in the slope of the power spectrum are far less im- 
portant in determining the formation history of a cluster 
than the actual value of ilo (Lacey & Cole 1993). Further- 
more, recent high-resolution simulations of the formation of 
dark matter halos in an Qo = 1 universe by Cole & Lacey 
(1995) demonstrate that a statistic similar to the one we use 
below to characterize cluster shapes is insensitive to spectral 
index and halo mass (cf their Figure 17) . Thus, our rescaling 
should have a negligible effect on our results. 

We tripled the size of our cluster sample by choosing 
three perpendicular axes through the centre of each cluster 
and treating the projection along each axis as a different 
cluster. This resulted in 72 clusters, 24 from each of the 
three cosmologies. 

2.4 CCD images 

To simulate typical observational datasets, we constructed 
artificial B-band CCD images of lensed field galaxies, build- 
ing them up pixel by pixel, in the manner described in WCF. 
We employed the same distributions of galaxy ellipticity, 
scalelength, redshift, magnitude, noise and seeing as in that 
paper: 

• Redshift and magnitude distributions 

For the redshift distribution of the model galaxies, we 
adopted the me = 25 distributio n predicted by th e ana- 
lytic model of galaxy formation of |Cole et al. (1994 ) . Since 
the critical density , Sa it, depends on the source redshifts 
through equation (2.2), we sampled this redshift distribu- 
tion discretely and produced a set of source planes spanning 
a range of redshifts. The net effect on the reconstruction 
equations is that Scrit is replaced by the mean value, Scrit, 
defined by 



s:. 



Scrit (z) 



p{z)dz, 



(2.10) 



where p{z) is the probability that a galaxy lies at red- 
shift z. For our adopted redshift distribution, Scrit — 6 x 
lO^^/iMo/Mpc^, with the exact value depending on the cos- 
mology. The distribution of apparent magnitudes we gen- 
erated directly from the B-band source counts of Metcalfe 



et al. (1995) 



Scalelength distribution 

We assumed that all the background galaxies are disks 
with exponential profiles. The scalelength was chosen from a 
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a) b) 

Figure 1. Examples of surface overdensity maps for four clusters from each of the three cosmologies. In each group the left-hand column 
shows the original surface density maps and the right-hand column the corresponding reconstruction. The greyscales have been adjusted 
according to the central surface overdensity in each case. From top to bottom the clusters are from the Hq = 1, open = 0.2 and fiat 
Qq = 0.2 models with the same initial fluctuation spectrum. The original CCD frame from which the reconstruction was made subtended 
an angle of ~ 10 arcmins but these maps show only the inner 6.7 by 6.7 arcmins. At the cluster redshift, z = 0.18, 1 arcmin is equivalent 
to 0.12/i~^Mpc for the Oq = 1 cosmology. The contour plotted encloses an area of 3.8 arcmin^. A Gaussian smoothing of width Ssm = 
0.25 arcminutes has been used in the KS reconstruction procedure. 



uniform distribution in the range 0.25 arcseco nds to 0.65 ar c- 
seconds, as suggested by the observations of Tyson (1994). 
• Elhpticity distribution 

Elhpticities for the background galaxies were generated 
by randomly sampling the empirical elhpticity distribution 
derived from a single frame in 0.7-0 .9 arcseconds seeing by 
Brainerd, Blandford, fc Small (199^ ). 

We created one CCD frame per cluster projection and 
then analysed it using FOCAS ( Jarvis fc Tyson 1981 ) as 
described in WCF. We assumed moderately good seeing of 
1 arcsecond full-width half-maximum. 



3 RESULTS 

Figures |^a-d show examples of the projected mass distribu- 
tion in the central regions of some of our clusters. In each 
group, from top to bottom, we illustrate clusters grown from 
identical initial conditions in the fio = 1, open Qq = 0.2, and 
flat Slo = 0.2 cosmologies respectively. The left-hand panels 
give the original surface density, smoothed with a Gaus- 
sian of Oam = 0.25 arcminutes, the same smoothing as we 
employed in the KS reconstruction. The right-hand panels 
show the corresponding reconstructions. The field of view 
in each plot is 6.7 arcmins (equivalent to 0.8h~^ Mpc, for 
ilo = 1), and the contour on each map encloses an area of 
3.8 arcmin^. 

The surface overdensities on the left-hand side have 
been normalised so that the mean value is zero over the full 
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Figure 1. continued. 



CCD frame, reflecting the fact that the mean value from 
the KS reconstruction is automatically set to zero. To en- 
hance the substructure in the reconstructed maps, we have 
used different greyscales in the left- and right-hand panels. 
The former appear uniformly black near the centre. Limita- 
tions of scale, particularly in the open ilo = 0.2 case, con- 
ceal the reality that the density distribution is in fact very 
highly peaked. The maximum surface density can be critical 
or even greater. In contrast, the maximum surface overden- 
sity in the right-hand panels is ~ 0.2 or 0.3 of critical. This 
diminution is partly due to observational effects such as see- 
ing and noise, and partly due to the effects of nonlinearity - 
weak lensing reconstruction techniques generally underesti- 
mate a in regions where the surface overdensity rises above 
a few tenths of critical. This failure of the method results 
in reconstructed profiles that saturate near the cluster cen- 
tre producing a broad plateau which, in the open flo — 0.2 
model, extends over many hundreds of pixels with similar 
density values. 

Figure 1 illustrates how clusters in the f2o = 1 cos- 
mology are still forming at z — 0.18. In most cases, large 
clumps are still infalling and the dominant mass conden- 



sation is usually quite elongated. By contrast, clusters in 
the open f^o = 0.2 cosmology are more centrally concen- 
trated and closer to spherical symmetry. Clusters from the 
flat Slo = 0.2 (Ao = 0.8) cosmology are similar to these 
but they have slightly lower central overdensity and are a 
little more elongated. In general, the reconstructed shapes 
match the originals rather well. With the exception of a few 
noisy pixels, the essential features of the different cluster 
morphologies are preserved by the reconstruction. 

We now us e th e dipole and quadrupole statistics, de- 
fined in Section 2.2, to quantify the distribution of cluster 



shapes in each cosmology. Initially we evaluated the statis- 
tics D{A) and Q{A) (equations 2.3 and 2.7) after identifying 
the cluster centre with the pixel of maximum surface over- 
density. We found that unless the reconstructions have very 
high signal-to-noise, this choice of centre results in noisy es- 
timates of D{A) and to a lesser extent of Q{A). For clusters 
with steep density profiles the reconstructed profiles have 
a saturated core and, as a result, the highest density peak 
is simply the highest noise peak anywhere in this central 
plateau. In light of this we tried a second choice of cluster 
centre which is more robust than simply selecting the dens- 
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Figure 2. The quadrupolc, Q{A), as a function of the area over which it is calculated. Two examples are shown of clusters in each of 
the Qq = 1 (left panel) and open Qq = 0.2 models (right panel). 



est pixel. This time we took the centre of coordinates to be 
the position at which D{A) = (i.e. the centroid of the area 
within the contour). The choice of centre now depends on 
the area, A, selected, but is much less sensitive to the pres- 
ence of noise. Having chosen the centre such that D{A) = 
we then use only Q{A) as a measure of the asymmetry of 
the cluster. 

Figure Ha shows Q{A) as a function of area. A, for two 
clusters from the fio = 1 cosmology. For A < 2 arcmin^, 
Q{A) tends to be small because of the dominant contribu- 
tion from the central relaxed core, although its value fluctu- 
ates due to noise. As larger areas (and hence smaller surface 
overdensities) are considered, some of the more distant in- 
falling clumps are included and Q{A) increases. Q{A) then 
varies slowly, gradually falling as an increasing number of 
randomly distributed noise peaks enter into the analysis. 
Fig. ^ shows Q{A) for two clusters from the open f2o = 0.2 
cosmology. Q{A) is fairly stable over a wide range of area, 
never rising above a value of 0.1. This behaviour reflects 
the smooth and nearly spherical mass distribution in these 
clusters. 

We now examine the distribution of Q{A), at a fixed 
area A, for our sample of clusters in each of the 3 cosmolo- 
gies. Figure 3 shows cumulative distributions of Q{A) ob- 
tained from the 24 clusters in each model (see Section |^. 
The three panels are for A = 3.8 arcmin^, 6.3 arcmin^, and 



12.6 arcmin . The two solid lines are for Q.n 



1, the two dot- 



ted lines for the open Oq = 0.2 model, and the two dashed 
lines for the flat fio = 0.2 cosmology. For each pair, the line 
marked by triangles, usually with the smaller quadrupole 
value, corresponds to the original, smoothed clusters. The 
unmarked line corresponds to the reconstructed clusters. 
Examining the original clusters flrst, it is apparent that an 
open 0,0 = 0.2 universe forms clusters with the smallest 
quadrupoles, followed by the flat 57o = 0.2 model and fi- 
nally the fio = 1 case. After reconstruction, the quadrupole 
value tends to increase, at least for the fio = 1 and the flat 
Slo = 0.2 clusters. This is primarily due to the inclusion of 
outlying noisy pixels in the analysis. 

Figure 3 shows that even with the imperfections of the 
weak lensing reconstructions, there remains a strong and 
measurable difference in the expected distribution of clus- 
ter shapes in high and low f2o cosmologies. The dependence 
of these distributions on Ao is very weak and so this statis- 
tic cannot be used to constrain the cosmological constant. 
This does mean, however, that f2o can be constrained in- 
dependently of the unknown value of Aq. For example, for 
A — 3.8 arcmin'^ we see that 75% of clusters in an f2o = 1 
universe have a Q[A) value in excess of 0.1, whereas less 
than 10% of clusters in the two fio = 0.2 universes are so 
aspherical. Alternatively, the median value of Q[A) is ap- 
proximately 0.2 for both A = 3.8 and 6.3 arcmin^ in the 
fio = 1 universe, but less than 0.06 for both the fio = 0.2 
universes. These large differences suggest that weak lensing 
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Figure 3. The cumulative probability distribution of the 
quadrupole, Q{A), at a fixed area, A, for A = 3.8 arcmin^ (top), 
6.3 arcmin^ (middle), and 12.6 arcmin-^ (bottom). The cluster 
centre has been chosen as the point about which the dipole van- 
ishes. The two solid lines are for the = 1 model, the two dotted 
lines for the open = 0.2 model, and the two dashed lines for 
the flat Ho = 0.2 model. In each case the curves marked by trian- 
gles correspond to the original clusters while the unmarked lines 
correspond to the reconstructed clusters. 



observations of a small number of clusters, approximately 
five or so, can distinguish between these two values of flo- 
A larger cluster sample is required to place finer constraints 
on SIq. 

We stress that the distributions of Q{A) plotted in Fig- 
ure 3 are typical of those expected from KS mass reconstruc- 
tions using weak lensing data obtained in realistic observing 
conditions. However, they are not intended to be definitive 
descriptions of Q{A) for these cosmologies. Clearly, the exact 
form of the distributions will depend on the observing con- 
ditions, the seeing and limiting magnitude, as well as on the 
redshifts of the clusters. In practice, simulations that mimic 
specific observational datasets will need to be performed for 
this test of fio to be reliably applied. 



4 DISCUSSION AND CONCLUSIONS 

We have shown how global properties of the dark matter 
distribution in clusters, as revealed by weak gravitational 
lensing analyses, can be used to measure the cosmologi- 
cal density parameter, f lp. Our approach com plements and 
extends earlier work by Evrard et al. (1994) who first ap- 
plied in practice the well established correlation between 
the morphology of clusters and the underlying cosmology to 
estimate the value of flo- Their analysis exploited the fact 
that the distribution of hot X-ray emitting gas in clusters 
is expected to refiect the structure of the underlying mass 
distribution. By probing this distribution directly, gravita- 
tional lensing bypasses the need to assume that the mor- 
phology of the gas faithfully mimics the morphology of the 
cluster mass. Possible (although unlikely) non-gravitational 
processes that may aff'ect this correspondence are therefore 
not a concern, nor is there a worry that the test may be 
affected by contaminating X-ray signals from, for example, 
AGNs. A second advantage of the lensing approach is that it 
allows the matter distribution to be probed at larger cluster 
radii than the centrally concentrated X-ray emission, that is, 
at radii where the distinction between different cosmological 
models is particularly strong. 

The major disadvantage of the lensing method is its 
sensitivity to projection effects. The observed galaxy shear 
pattern is a function of the product of the projected sur- 
face density along the line-of-sight and the effective cross- 
secti on f or lensing (which depends on distance through equa- 
tion 2.2). Mass clumps at large distances from the cluster 
contribute to the lensing signal and, in principle, to the mea- 
sured quadrupole. In practice, this is unlikely to be a strong 
effect, particularly when the analysis is restricted to rela- 
tively small areas around the cluster centre. Nevertheless 
the size of this contamination needs to be assessed using 
different simulations from those discussed here. X-ray anal- 
yses are much less sensitive to projection effects and so the 
two approaches are complementary and should be used in 
combination. 

In some respects, our detailed r esults must be regarded 
as preliminary. The simulations by Evrard et al. which we 
have analyzed have a number of limitations. Firstly, the ini- 
tial conditions in all the cosmological models were laid down 
with the power spectrum appropriate to an fio = 1 CDM 
universe. This inconsistency between the power spectrum 
and the cosmology was further exacerbated in our case by 
the need to rescale the simulations in order to ensure a lens- 
ing signal with comparable signal-to-noise in all cosmolo- 
gies. Secondly, the baryon fraction in all the simulations 
was fixed to be 10% of the critical density so that in the 
flo = 0.2 models, the gas contributes half of the total grav- 
itating mass. Finally, the clusters we have analyzed were 
chosen to correspond to high peaks in the density field, but 
they do not represent a proper statistical sample. Although 
we argue that all these approximations are unlikely to affect 
our results significantly, it is clearly desirable to repeat our 
analysis with purpose-built simulations. 

The main result of our work is that the distribution of a 
low-order statistic sensitive to global deviations from spher- 
ical symmetry in the mass distribution of clusters depends 



8 G. Wilson, S. Cole, and C. S. Frenk 



strongly on Q,o . The distribution of this quadrupole statistic 
can be robustly recovered from the shear pattern of galax- 
ies weakly lensed by the cluster gravitational potential. We 
have shown that under standard observing conditions, lens- 
ing data for a handful of clusters should distinguish between 
cosmological models with f2o = 1 and fio = 0.2. 
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